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Background 

Different methods have been proposed to estimate Wall 
Shear Stress (WSS). Morgan and co-workers (J Thorac 
Cardiovasc Surg and Ann Biomed Eng 1998) used a 
finite difference scheme to quantify the WSS tensor. 
However, it is well known that finite-difference methods 
cannot effectively handle complex geometries, as those 
found in the cardiovascular system. To account for arbi- 
trary cross-section shapes, Stalder et al., MRM 2008, 
used B-spline (BS) interpolations to smoothly describe the 
lumen contours. In this work, we propose and validate a 
new method for calculating the WSS distribution based on 
Finite-Element (FE) interpolations. 

Methods 

The velocity field obtained at discrete locations from 
2D and 3D CINE PC-MRI was interpolated using lin- 
ear triangular FE. The shear stress tensor over the 
entire section cut was obtained from a global least- 
squares stress-projection method, from which the axial 
WSS vector was obtained. The proposed approach was 
benchmarked against a modified Poiseuille flow profile, 
and the robustness of the method was assessed by 
changing the level of resolution and noise. Addition- 
ally, we computed the WSS distribution in different 
aortic sections from a pulsatile aortic phantom, and 
from 5 healthy volunteers. In the aortic phantom and 
volunteers, 3D CINE PC-MRI flow data was acquired 



in 2D cutting planes in 5 with a spatial resolution 0.8 
mm and temporal resolution 37 ms. We have also 
compared our framework with a BS based method 
previously reported in the literature. 

Results 

Our results showed that the local WSS values were in 
good agreement with the theoretical values obtained 
from the modified Poiseuille flow problem. The aver- 
aged WSS over the vessel contour showed a systematic, 
but negligible bias compared to the Poiseuille averaged 
WSS when subjected to different levels of noise and 
resolution, (see Figure 1). In contrast, the BS based 
method led to greater local differences, and in averaged 
WSS values that were largely affected by the level of 
noise and resolution. In volunteers, the cardiac cycle 
average value of the average WSS value was 0.21 ± 
0.06 N/m^, whereas the BS based method yielded 0.45 ± 
0.13 N/m^. The Bland-Altman plot (Fi gure 2) showed a 
systematic bias between both methods with a mean 
WSS difference of -0.2338 ± 0.2491 N/m^. 

Conclusions 

In conclusion, we have developed a novel methodology 
to calculate WSS based on FE interpolations, which 
provides an excellent approximation of local WSS 
values, stability when subjected to noise and remarkable 
convergence properties as the pixel size is decreased. 
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Figure 1 WSS contour mean and standard deviation (left column), and root mean square error (RMSE) (right column) for a noise level 
of 0% of maximum velocity peak (A - D), 0.05% (B-E), and 0.1% (C - F). The methods analyzed were the modified Poiseuille flow (P), 
the proposed method (FE), and the B-Spline based method (BS). 
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Figure 2 Bland-Altman plot of cardiac phase averaged wall 
shear stress contour mean comparing the proposed method 
(FE) and the B-spline method (BS) from volunteer data. 
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